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Abstract 

This work is devoted to the investigation of possibility of controlling of ions motion inside Paul 
trap. It has been shown that by proper selection of the parameters of controlling electric fields, 
stable localization of ions inside Paul trap is possible. Quantum consideration of this problem is 
reduced to the investigation of the Mathieu-Schro dinger equation. It has been shown that quantum 
consideration is appreciably different from classical one that leads to stronger limitations of the 
values of the parameters of stable motion. Connection between the problem under study and the 
possibility of experimental observation of quantum chaos has been shown. 

PACS numbers: 73.23.-b,78.67.-n,72.15.Lh,42.65.Rc 
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I. INTRODUCTION 



It is well known that formation o 



miniature trap technology 



' Bose-Einstein condensate is based on application of 



Formation of condensates of alkaline atoms turned out to be possible thanks to unification 
of new ideas and technologies developed in other areas of physics. First of all magnetic trap 
must be formed in order to confine the condensate. 

In spite of this, the traps formed by application of variable electric fields are the subjects 
of special interests. The point is that, the confinement of separate ions in trap for a long 
time gives new possibilities in laser spectroscopy. In addition, single entrapped ion represents 
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itself a unique object to prove the fundamental laws of quantum mechanics 13, |8| 
For example, single ion placed inside Paul trap is applied to realize a quantum gate, and a 
chain consisting of many ions pla ced in linear trap may be regarded as promising component 



of future quantum computer 



12| 



T he g oal of this work is the study of quantum dynamics based on the application of Paul 



trap [l3| 



OH, 



16] 



Not going into details that can be found in [16j, we shall describe shortly necessary 
information about Paul trap needed further. 



Paul trap is based exceptionally on the application of sign-varia 
potential formed by electrodes inside trap has the following form 



Die electric field. Besides, 



2^2 



(r,t) = f(t)(x 2 + y z z 



2z< 



(1) 



where 



f(t) 



U + V cos(u r ft) 
rl + 2z 2 



Here, U and V are the amplitudes of constant and variable fields, ro-radius of circular 
electrode, 2z -distance between covers, and u r f = 2ir/T frequency of variable field. 
The ion's classical motion in potential (1) is described by: 



Mr — -eV0(f,t), 



(2) 



where M and e are mass and charge of the ion. Putting potential (1) into equation of motion 
(2) yields: 

M?+2ef(t)(r-3z) = 0, (3) 



where z = (0, 0, z). 
If we label: 



ip — x, If 



1 



8eU 



2^' a ~ Mul f {rl + 2ziy 



AeV 



2ziy 



(4) 



Eq. (3) may be rewritten in the form of the Mathieu equation 16[ : 



dtp 



+ (a + 2q cos(2(p))ijj(ip) = 0. 



The Mathieu equation is a well known equation of mathematical physics 



usually encountered when studding the problems of parametric resonance 18 



(5) 

and it is 



The presence of the areas of stable and unstable motion is a peculiarity of parametric 
resonance. The solutions of Eq. (5) may be as stable as unstable depending on the values of 
a and q. In case of Paul trap, ion confinement inside a trap corresponds to stable solution. 
Because of this, we are interested in the stable solutions (2) of Eq. (5). Depending on 
physical situation for studying motion of ions inside Paul trap classical and quantum con- 
siderations are used. The modern state of affairs in studying this problem is given in 16] . 
In spite of the fact that, classical consideration of this problem is more or less adequate, 
quantum-mechanical consideration is far from perfection 16) . 

The goal of this work is to fill up this gap. 

In this work we shall try to give more adequate (in our opinion) quantum mechanical 



description of motion of ions inside Paul trap. At the same time we shal 



obtained by us from the study of the problem of quantum chaos 
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use the results 



23|]. 



II. STABLE QUANTUM-MECHANICAL SOLUTIONS. THE MATHIEU- 
SCHRODINGER EQUATION SYMMETRIES 

The Mathieu Eq.(5) in classical consideration , depending on the values of the parameters 
a and q, has as stable as unstable solutions. On the plane (a,q), the boundaries between 
stable and unstable solutions pass through special curve called the Mathieu characteristics 
[17]. Along these lines, to which certain set of the values of the parameters (a(q),q) cor- 
respond, the well known periodic Mathieu functions ce n (tp,q), se n (ip,q) are solutions 
of classical Mathieu Eq.(5). Periodic solutions are the subject of our special interest, since 
periodic solutions corresponds to finite motion, which means confinement of ion inside Paul 
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trap. It is worth noting that, this kind of selection of the values of the parameters of a(q) 
and q impose limitations on the amplitudes of constant and variable fields: U and V (see 
Eq.(4)). 

In quantum consideration situation gets considerably complicated. The point is that, in 
quantum case Eq. (5) must be considered as the Mathieu-Schrodinger equation. In this 
case solution ip(tp,q) and the value of parameter a(q) must be considered as eigenfunction 
and eigenvalue of quantum-mechanical equation. Characteristic property of the Mathieu- 
Schrodinger equation is the following: eigenfunctions ip n (<p, q) can differ from the solutions 
of classical Mathieu equation ce n (<p,q), se n (<p,q). Finding of eigenfunctions ^ n (y?,g) of the 
Mathieu-Schrodinger equation is difficult problem and the application of group theory is 



2l|. In the present work we shall only 



needed. Interested reader can find details in work 
use the results needed further. 

The eigenfunctions of the Mathieu-Schrodinger equation form the basis of irreducible 
representation of the Mathieu-Schrodinger equation symmetry group. 

As is known [17], periodic solutions of Eq. (5) are given by the Mathieu functions: 

ce 2n (f,q), ce 2 n+i(y?,g), se 2n +i(f, q), se 2n +2(f,q) (6) 
which setisfy the normalization condition: 
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where ijj n ((p,q) means Mathieu functions (6). To eigenfunctions (6) there corresponds the 
eigenvalues (Mathieu characteristics) 

(?), (?), (7) 

which also depend on the parameters q. 

The properties of symmetry of the Mathieu functions can be presented in the form of 



Table (see Table 1 in Ref. |2l|). 

By immediate check it is easy to be convinced that symmetry properties of the Mathieu 
functions with respect to four elements of transformation: 

G(<p — > —ip) = a, G(ip — > it — (/?)= b, 

G{ip — >• 7r + <f) = c, G((p — > tp) = e 



form a group. For this purpose it is enough to test the realization of the following relations: 

2 l2 2 

a = b = c = e, 



ab = c, ac = b, be = a. 

Group G contains three elements a, b, c of the second order and unity element e. The 
group G is isomorphic to the well-known group of Klein 24j. 

The group of transformations G is not a simple group since it contains subgroups. When 
combined with the unit element, each of three elements a, b, c forms a subgroup of second 
order: 

G + — > e, b; 
G_ — > e, c; 

Go — ► e, a. (8) 

Thus the Mathieu-Schrodinger equation is characterized by a presence of three symmetries 
determined by subgroups of symmetry group. Depending on the values of parameter q, the 
system can be in the states characterized by any of the symmetries (8). It is worth mention 
that, to the symmetries G + and G_ correspond degenerate states with wave functions [21]: 



a/2 

G - ^fn+M = ^-(ce 2 n+i(<^) ± ise 2n+ i(<p)), 



V2 

Tp2n(<P) = ^-(ce 2 n(<^) ± ise 2n (ip)) (9) 



and 



a/2 

G+ -> V2n(v) = ~y ( ce 2n(<P) ±ise 2n+1 (<p)), 
a/2 

vL+M = ^-(ce 2 n+i(<^) ± ise 2n +2{^))- (10) 



To the symmetry group G$ correspond nondegenerate states with wave functions 



2l|: 



G -> ce 2n { l f); ce 2n+1 (ip); se 2n ((f); se 2n+1 ((p). (11) 
In the end of this section, let us sum up some results. 

The Mathieu-Schrodinger equation is characterized by a specific dependence of the spec- 
trum of eigenvalues a n (q) and eigenfunctions ip n (tp,q) on the parameter q (see Fig.l). On 
the plane with the spectral characteristics (so-called Mathieu characteristics [17]) of the 



problem, this specific feature manifests itself in the alternation of areas of degenerate (G±) 
and no degenerate (Go) states. The boundaries between these areas pass through the branch 
points of energy terms a n (q). 

In case of quantum consideration, the problem of confinement of ion inside Paul trap, the 
aforesaid means the following: 

In order ion's position to be localized inside Paul trap, it is necessary that the solutions 
of the Mathieu-Schrodinger quantum-mechanical equation to be periodic and as a result 
finite. This can be achieved only by proper selection of the values of the parameters a(q), 
q corresponding to the Mathieu characteristics, since along these lines motion is periodic. 
According to Eq. (4), this can be done by proper selection of the parameters of electric field 
U and V. The system with these kind of selected values of the parameters can be found in 
the state characterized by one of the symmetry and it shall be described by corresponding 
wave functions (9)-(ll) (see Fig.l). 

In case of classical consideration the situation is different. Not only the Mathieu charac- 
teristics belong to the areas of stable motion but also the areas where the Floquet solutions 



are damped functions of time 25]. 



In case of quantum consideration the areas of damped solutions belong to forbidden 



energy zones in semiconductors 26[. This is characteristic property of quantum case. 



III. NONSTATIONARY CASE. STOCHASTIC HEATING OF THE SYSTEM 

Let us assume, that the amplitude of the electromagnetic field V is modulated by a slowly 
changing field. According to (4), the influence of modulation can be taken into account by 
making a replacement in the Mathieu-Schrodinger equation (5) 

q(i) = <?o + Agcosz/t, (12) 

where Aq is the modulation amplitude expressed in dimensionless units and v is the modu- 
lation frequency. In this case, the problem of ion's motion inside Paul trap is reduced to the 



problem studied in work 22|. In the mentioned work stochastic absorption of energy by non 



stationary chaotic quantum system was studied. Not going into details that can be found 



by interested reader in 22| , here we only present the main results related to given problem. 



After making replacement (12), Hamiltonian for the Mathieu-Schrodinger equation takes 
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the form: 

H = H + H'(t), 
d 2 

H = —q-^ + <?ocos2</?, 

H'(t) = Agcos2<pcosz/i. (13) 

We assume that gradual change of q(t) may involve some N branch points on the left and 
right side of the separatrix a(q) = q, (see Fig.l). 

Simple calculations show that the matrix elements of perturbation H'(t) with respect to 
the wave functions of the nondegenerate area Go (11) are equal to zero 

< ce n \H'(t)\se n >= 0. (14) 

But in the degenerate areas G + , G_ 

H' + _ =< ^ 2n+1 \H'{t)\^ n+l >^ 0, (15) 

where ip 2n+1 axe the functions (9), (10) corresponding to the areas of symmetry G_ and G + 
(see Fig.l). 

Thus the modulation (12) may cause transition between quantum states of the degenerate 
areas G + , GL. 

When there is a small dispersion in the values of parameter 5qo , as was shown in work 



22l | . mixed states are formed in the degenerate areas. As a result of this, the system 
may be with equal probability of 1/2 in two states ip2n an d i>2n- Consequently, because 
of specific dependence of energy terms a n (q) on the parameters q (see Fig.l), and because 
of time modulation q(t), the process is extended to other levels as well. At the same time 
absorption of energy by the system and stochastic heating takes place (see Fig. 8 of work 



Something like that must take place in case of ion placed in Paul trap and kept there by 
amplitude modulated alternate electric field. 

Thus Paul trap may be used for experimental observation of quantum chaos. The present 
theoretical analysis gives some recommendations concerning the values of dimensionless 
parameters of the Paul trap(4): a pa 6.7 -j- 17.2; q pa 12. 

The point is that, the problem of internal rotational motion in polyatomic molecules is 



reduced to the Mathieu-Schrodinger equation model 



271 ] . Therefore as is shown in work 



FIG. 1: A fragment of the parameter-dependent energy spectrum. Squares denote branch points. 



23j, the present experimental data of infrared absorption in molecules CiH§ 



are direct 



evidence of proposed model to be correct. 
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